Current-Induced Motion of Narrow Domain Walls and Dissipation in Ferromagnetic 
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Spin transport equations in a non- homogeneous ferromagnet are derived in the limit where the 
sd exchange coupling between the electrons in the conduction band and those in the d band is 
dominant. It is shown that spin diffusion in ferromagnets assumes a tensor form. The diagonal 
terms are renormalized with respect to that in normal metals and enhances the dissipation in 
the magnetic system while the off-diagonal terms renormalize the precessional frequency of the 
conduction electrons and enhances the non-adiabatic spin torque. To demonstrate the new physics 
in our theory, we show that self-consistent solutions of the spin diffusion equations and the Landau- 
Lifshitz equations in the presence of a current lead to a an increase in the terminal velocity of a 
domain wall which becomes strongly dependent on its width. We also provide a simplified equation 
that predicts damping due to the conduction electrons. 



Dynamics of magnetic domain walls (DW) is a classic topic [l], 0, 0] that recently received a lot of attention 
i— | ■ due to new fabrication and characterization techniques that permit their study at the nanometer scale. Moreover, 
the subject of spin dynamics in the presence of large inhomogeneities is currently of great interest experimentally 
i-^h ■ and theoretically due to the potential applications in various nano-devices, especially magnetic storage 4|. One 
particular area that is still not well understood is the interaction of domain walls (DWs) with polarized currents. 
■ The question here is how best to represent the contribution of the spin torque to the dynamics of the magnetization 
SBSBBES II, 12) 13, 14, Il| • So far attention has been focused on wide DWs where it was shown that terminal 



velocities are independent of the DW width [3, [lj] . 

This paper extends previous treatments to the case of thin, less than 100 nm, DWs. One of the main objectives 
of our work is to expose the interplay between linear momentum relaxation and spin relaxation as the conduction 
electrons traverse a thin DW. This interplay originates from the strong exchange interaction between the conduction 
s electrons and the localized d moments, and makes the terminal velocities as well as the transport parameters of 
the conduction electrons dependent on the configuration of the local magnetization. This leads to an enhancement 
: of the non-adiabatic contribution of the spin torque to the DW motion and opens the way to study spin torque- 

induced magnetization dynamics in thin DWs in greater depth by measurement of DW velocities. Moreover, we show 
t— I ■ that the interaction of the conduction electrons and the d moments is also relevant for homogeneously magnetized 
metallic systems, where it is at the origin of intrinsic damping. Our work can be easily adapted to magnetic multilayer 
structures and hence the equations derived here are capable to treat non-collinear magnetization geometries as opposed 
to that in ref. flil w hich deal only with collinear configurations. Narrow DWs can exist either naturally [It], EH 
or artificially 13 ? 20| and we hope the results discussed here show the potential benefits of studying dissipation 
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DW-like structures. 

To derive the spin coupling of the s electrons to the magnetization, we adopt the sd picture which has been the 
basis for most of the studies in DW motion In the following we use (I, m, n) for moment indexes, and k) for 
■ space indexes. In addition the transverse domain wall is assumed to extend in the x direction, with magnetization 
\ pointing in the z direction. We start from the Boltzmann equation satisfied by the 2x2 distribution function of the 
conduction electrons, f = f e + f s • a, where er; (I — 1, 2, 3) are Pauli matrices, in the presence of the magnetization M 
of the system and an external electric field E: 



d t { + v ■ Vf +e (E + v x H) • V p f + 
i[ M sa-H sd ,f] =- J - *° °. (1) 

Tp T sf 

The sd exchange field is H s rf(x, t) = JM(x, t)j [is with J sw 1.0 eV, and t p , T s f are the momentum and spin relaxation 
times, respectively [lil l2ll. |22|. The variables v, e, and [Ib are the velocity, the charge and the magnetic moment of 
the s electrons, respectively, f® and f° are the equilibrium charge and spin distribution. 

The conduction electrons have a polarization m = [Ib J j^jTraf and carry a charge current j c = e J j^jvTrt, 
as well as a spin current 

^vTraf. (2) 
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In the following we use normalized definitions of the moments, i.e. 
assumed to satisfy a Landau-Lifshitz-Gilbert (LLG) equation 



dM 

dt 
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1. The d electrons will be 



(3) 



where r ex is the inverse of the precessional frequency, uj c — J/h, of the conduction electrons due to the exchange 
field. H e ff is the total field acting on the magnetization which includes the exchange field between the d- moments, 
the demagnetization field and the anisotropy field. In metals, the main source of dissipation is believed to be due to 
the conduction electrons which in our theory is accounted for explicitly within the limitations of the sd model [23| . 
Hence, the damping constant a p d is assumed to be due to dissipation caused by channels other than the conduction 
electrons such as phonons or defects. 

In inhomogeneous magnetic media, the sd exchange term becomes comparable to that of the Weiss molecular field 
and hence the effect of the conduction electrons on the magnetization should be taken beyond the linear response 
approach. Going beyond the linear theory will allow us to see how the presence of the background magnetization 
affects the transport properties of the conduction electrons. We believe this is especially true in transition metal nano- 
magnetic devices where the hybridization of the s and d electrons is strong. Using standard many-body methods [l5| , 
the diffusion contribution to the spin current can be found 



j l : (t,x) = -P to (t,x)V i m n (t,x), 



(4) 



where T> is a diffusion tensor with effective relaxation time r which will be assumed equal to the momentum relaxation 
(t ps t p ). The V tensor obeys the reduced symmetry of the ferromagnetic state and is (lBj 
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where Vt = JM/h, D A 



with D 



2 t p being the diffusion constant of the electron gas with Fermi 



velocity Vf. It should be observed that in the presence of spin-orbit coupling, the symmetry of the diffusion tensor 
will be the same as given here but the separation of the relaxation times in independent channels of momentum and 
spin relaxation will not be valid. In the following, the effect of the electric field is taken only to first order. 

The symmetry of the spin current is best revealed by going to a local frame where the magnetization lies in the 
z-direction. In this frame, one obtains for E = 



3± 



dm 
dx ' 
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where m (x) 



(x) 



i y (x), and D c ff = D± + iD xy is an effective diffusion coefficient with D xy = D±tlo c . From 



the divergence of the spin current we get the steady-state equation for the spin accumulation, 



dx 2 
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where Ag ff = T e gD c s with t c s = l/(^- — ioJ c ), m is the equilibrium spin density, and X s( u is the longitudinal spin 
diffusion length typically in the range of 5-100 nm. The general solutions for the complex accumulation are of the 
form m{x) = Aexp [— x/A e ff] + -Bexp[a;/A ff], i.e. they show an exponential decrease (or increase) and oscillations 
from a local inhomogeneity in M. In the limit of a large sd exchange field the period of the oscillations is ^ which 
corresponds to the coherence length l/\k' — k^\ in the ballistic approach, where is the spin- up momentum. 

Our expressions for the spin current generalize those used currently in the literature |9| . We find that the diffusion 
constant Do is now renormalized by 1/(1 + (tu c ) 2 ) which means that precession in the exchange field reduces diffusion. 
Moreover, the precession gives rise to off-diagonal terms in the diffusion tensor which reflect the local 2D rotational 
symmetry around M. 

The origin of the off-diagonal term T> xy can be understood qualitatively in terms of flux. First we rewrite it in the 
following form 



D xy = (-v f t p )- 



Vf 2 W c 



iv 2 



(8) 
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where v = In the limit of fast precession, v <§C co c , we have D xy = ^V[ 2 /ui c . Next if we set V[/ui c = L m , then L m 
is the distance a spin typically goes before it 'converts' into the spin at 90-degrees to that which it started with. The 
corresponding contribution to the flux has an obvious interpretation - the source of spin x, m x , is particles coming 
from a distance L m away where they had spin y, m y . The flux can be derived from a simple 'kinetic' argument. A 
distance L m upstream, the density is m y = m y + and a distance L m downstream, m y — m? — L„,^. The 

flux of particles with spin x, m x , crossing a point, coming from upstream, is m y ^Vf/3 and from downstream it is 
m y ^V{/3. The difference is then (m y ^ — m y ^)vf/3 = 2L m ^^Vf/3 which, within a factor of 2, is our off-diagonal flux. 
In short, the off-diagonal terms are the corrections induced by precession on the diffusion process. 

To get the effective equation for M, we use equation 3 to express m in terms of the magnetization M, then use 
the implicit solution back into the equation for m, Eq. 1. We find that the equation of motion for the magnetization 
becomes 

(3^ = - 7 M x H e// + a • VM + (a pd -px^ (9) 
-(7M x (M x H e //) + V • PVm, 

where a = P(^s/e)Trj, /3 = 1 + m + a p d£, (or MB/3' = (1 + m + £ 2 ) ), £ = T ex /r s f is the ratio of the precessional 
time to the spin relaxation time of the conduction electrons, and P is the spin current polarization Q. The second 
term on the right is the adiabatic spin torque while the last term is the diffusion contribution. The third and the 
fourth terms are equivalent to the non-adiabatic spin torque with the original damping a p d from Eq. 3 added to the 
third term. For uniform magnetization, VM = Vm = 0, and damping constant a p d = 0, Eq. [9] reduces to 

dM C 2 + /? A/r „ 

~df = -^h^ m x ff 

-7Z-0STp M x ( M x Hcff ) • ( 10 ) 

Hence, we are able to predict damping due to conduction electrons, quantify the corresponding damping constant 
«ei = 7C TO o/(/3 2 + £ 2 ), and identify its origin as the spin torque contribution of the conduction electrons. We have 
written Eq. [TO] in the LL form, but it equally can be written in the LLG form. We have already shown in ref.[l5| that 
the magnetization dynamics of a thin film embedded between two normal conductors and subjected to an electric 
field is not well described by closed LL (or LLG) equations. 

Next we discuss qualitatively the effect of the diagonal and off-diagonal terms of the diffusion tensor on the velocity 
of a domain wall, of width A. If we ignore the spatial dependence of the diffusion tensor elements and replace the 
Laplacian in the diffusion equation by 1/A 2 , then we recover equations similar to those discussed by Zhang and Li 
Q but with renormalized spin flip scattering rate, 1/t s / — * I/t^j = 1/t s / + Dq/X 2 , and renormalized precessional 
frequency, l/r ex — > = l/r ex — D xy /X 2 . Therefore, the velocity and the effective damping of the DW dependent 

on the size of the inhomogeneities in the magnetization. This can be understood qualitatively from the results in Q 
which showed that the DW velocity v for a wide DW, i.e. Vm 0, is inversely proportional to the damping a e i (in the 
case a p d—0), v (Pjfi B /e)((l + £ 2 )/(£ m o))- Then, ignoring the renormalization of the diffusion coefficient D n , the 
velocity is expected to take a similar form as in the case which does not account for the diffusion but with £ replaced 
by £ N — T^/rjj. The damping a will be also affected by this renormalization as is expected, since broadening due 
to inhomogeneities is well known to occur in ferromagnetic resonance measurements. 

Now, we turn to the discussion of the results of the above theory for a 1-D DW configuration. We solve numerically 
the coupled equations of motion for the conduction electrons and that of the magnetization. We include the d-d 
exchange between the local moments, the anisotropy along the direction of the current and the dipole field. Pinning is 
neglected but can be easily included in the simulations. Besides varying the width of the DW, we also vary the other 
parameters in the sd model since there is no universal agreement on their exact values. For example, it is generally 
believed that spin relaxation times are about two orders of magnitude longer than momentum relaxation times. While 
this may be true in paramagnets , we already know that in NisoFe2o they are comparable [24] . In Permalloy, the spin 
diffusion length, l s — Vf^/r s fr p , is of the order of 5 nm which is of the same order as the mean free path, l p = 3 nm. 

Figure 1 shows the effect of introducing the (unnormalizcd) diffusion term Dq in the equations of motion of the 
magnetization. For DW width larger than 100 nm our result approximately recovers that of Ref. 0]. The variations 
of the domain wall velocity v with A are found to depend strongly on Dq. This is expected since v is, to first order, a 
function of Dq/X 2 (cf. inset). Moreover, the velocity peaks when the mean free path of the conduction electrons, l p , 
is of the same order as the DW width, since for l p » A there is almost no scattering while for ! p «A there is only 
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FIG. 1: Domain wall velocity as a function of domain wall width for j c = 10 8 A/cm 2 , r s f = 3.0 x I0~ 12 s, a p d = 0.01 and 
three different diffusion coefficients, T> = DqI, given in the figure in units of m 2 /s. A value of Do = I0~ 2 m 2 /s corresponds to 
t p ~ I0~ 14 s. The solid thick line is that of Zhang and Li 0. The inset shows the DW velocity versus 

slow diffusion. We have confined our results to A > 2 nm since at much smaller DW widths, we expect contributions 
from Coulomb interactions and a breaking of the quasiclassical picture employed here. 

In figure 2, we show the effect of the corrections introduced by the off-diagonal terms in the diffusion tensor. This 
non-adiabatic effect actually appears to suppress the DW velocity or the effect of diffusion as we explained earlier. 
Otherwise, the functional behavior of the velocity remains similar to the one discussed in figure 1. 



t — 1 — i — ■ — i — ■ — i — i — i — 1 — r 
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FIG. 2: Domain wall velocity as a function of domain wall width with the correct diffusion tensor taken into account. The solid 
(open) symbols are without (with) off-diagonal corrections of the diffusion tensor. Parameters are identical to those in Fig.Q] 



Finally in figure 3, we extract the contribution of the conduction electrons to the effective damping of the magneti- 
zation. First, we observe that the off-diagonal diffusion terms have little effect on the relaxation of M which is mainly 
determined by the spin relaxation time r s f. These results are also not sensitive to the DW width and the extracted 
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electronic damping has the correct order of magnitude for metals. 
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FIG. 3: The electronic damping a e i as a function of spin flip scattering r s / for a 10 nm domain wall. The solid (open) symbols 
are for off-diagonal terms included (not included). The diffusion constant is Do = I0 _2 m 2 /s and a p d = 0. 

In summary, we have solved the conduction electron-magnetization problem in the presence of a current self- 
consistently. We found that the diffusion term provides a larger contribution to the drive torque than to the damping 
process, leading to a overall increase of domain wall velocity. We also showed that the new off-diagonal terms of 
the diffusion tensor enhance the DW velocities which become at least one order of magnitude larger than previously 
found. Moreover, the dependence of the DW velocity on the width of the DW was found to be non-linear and strongly 
dependent on the non-adiabatic behavior of the conduction electrons through the non-diagonal corrections of the 
diffusion tensor. We have been also able to determine the contribution of the conduction electrons to the damping 
in ferromagnetic metals which we found to be of the same order as the typical measured values of a. Therefore, 
our treatment allows us to include electronic damping in micromagnetic calculations in a more rigorous way than is 
currently done by simply accounting for it by a simple a parameter. 
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